Centre de Physique Theorique - CNRS - Luminy, Case 907 

F- 13288 Marseille Cedex 9 - France 

Unite Propre de Recherche 7061 

^ A STOCHASTIC PROCESS 

a^ : FOR THE DYNAMICS 

OF THE TURBULENT CASCADE 



(N 



in 



o 






c^ 



R. LIMA and R. VILELA MENDESQ 
Abstract 



Velocity increments over a distance r and turbulent energy dissipation 

on a box of size r are well described by the multifractal models of fully 

developed turbulence. These quantities and models however, do not involve 

O ' time-correlations and therefore are not a detailed test of the dynamics of the 

^ ■ turbulent cascade. 

If the time development of the turbulent cascade, in the inertial range, 

^! is related to the lifetime of the eddies at different lenght scales, the time 

'^ I correlations may be described by a stochastic process on a tree with jumping 

O I kernels which are a function of the ultrametric (tree) distance. We obtain 

the solutions of the Chapman- Kolmogorov equation for such a stochastic 

O ■ process, with jumping kernels depending on the ultrametric distance, but 

with an arbitrarily specified invariant probability measure. We then show 

k>< I how to use these solutions to compute the time correlations in the turbulent 

H ' cascade. 
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1 Introduction 

One of the most interesting phenomena in fully developed turbulence is the 
occurence of an energy cascade from the macroscopic length scale L of the 
experimental apparatus down to smaller and smaller length scales. Lenght 
scales / in the range L » I >> rj, rj being the scale where the fragmentation 
process is stopped by dissipation, are said to be in the "inertial range". 
In the inertial range viscosity effects are not important and Kolmogorov^ 
proposed long ago a scaling theory with conserved energy transfer between 
length scales. From scale invariance and the assumption that turbulence is 
space-filling it follows that the velocity fluctuation 6v{l) over an active eddy 
of size / scales as 

{\\sv(i)r) - 1^^ (1.1) 

with c = f . 

However turbulence may or may not be space-filling^ and the volume of the 
active eddies may change when the energy is transferred from the scale In 
to the scale In+i- This leads naturally to a fractal structure for the cascade 
with fractal dimension less than 3. For example in the /9-model^ the rate of 
energy transfer 



E, 



n 



'"n 

does not change along the cascade but the total mass of the active eddies is 
multiplied by (3 at each step. Then the exponent Cp in Eq-( 1-1) becomes 



Cp = p./i + 3 - D 
with h = ^^^1 where D is the fractal dimension related to (3 by 

log2/? = D-3 



)-n 



if the length scales are related by /„ = Iq2' 
The /5-model as well as a log-normal model^ for the distribution of En are 
however in contradiction with the experimental results on moments of higher 
order for the velocity structure functions^'^. This fact led to the proposal of a 
multifractal generalization^'^ called the random /?-model where it is assumed 
that, at each scale /„, there are several distinct /3„(/c)'s which are chosen 



according to some probability law. That is, the energy transfer may take 
place according to several distinct dimensional routes. Requiring a fixed 
energy transfer rate one obtains 

^"^^^^-/5.,,(fc)M±iM (1.2) 






ra+1 



Hence at scale /„ the velocity fiuctuation in each eddy depends on the frag- 
mentation history which is defined by the product Pi./32-'"-(3n- Then 



svn--a{Y[(3^r' (1-3) 

1=1 

and 

{\\SVn{Q\n-£jf[dMl''P{Pu-\Pn) (1.4) 

P(/5i,'" , /?„) being the occurence probability of the sequence j3i;" , /3„. 

At the level of precision of the existing experiments, agreement with the data 

is already obtained if one assumes independent fragmentations 

n 

P{(3l;■\Pn) = l[P((^^) (l-5a) 

and a simple binomial process 

F(/3) = 7(5i + (1 - 7)(5i (1.5b) 

7 is a parameter chosen to fit the data (7 ^ 0.875). 

The assumptions (1.5) only define the probability distribution pi at each 
level of the cascade tree (Fig.l). They make no statement concerning the 
time evolution and the time scales of the eddies in the cascade. In fact it 
is well known^'^ that only a few restrictions are imposed on the form of the 
velocity fields by the predictions of the statistical models described above. 
These are essentially the existence of singularities in the derivatives of the 
velocity field at some points. Aside from that, several distinct velocity fields 
may be compatible with the spectra and the scaling laws. They range from 
the superposition of random uncorrelated Gaussian components having only 
a correct spectrum in the inertial range^°'^^ to simple fiow fields in isolation, 



such as a vortex sheet wrapping up while being stretched by a large-scale 
straining motion^^. 

One way to further the research on the dynamical properties of the velocity 
field is to study time-correlations for the observables. To make a connection 
of the models with the time dependence of each observable one should first 
consider the dynamical aspects of the energy cascade itself. To do this a time 
hierarchy in the development of the turbulent cascade must be defined. In 
what follows we deal with this issue. 

For a binary cascade tree (Fig.l) we may use a dyadic labelling for the possi- 
ble states at each level. The state space Vn at level /„ is the set of all products 
^(i)_..._^(n) ^i^i^ p(i) ^ ^^^^^_ (jj^ Eq.(1.5b) po = 1 and pi = |). There are 

2" elements in Vn and the probability of the state i is 

p. = P{j3^^\- .[5^'')) = 7«o«(l - 7)"i« (1.6) 

where no(i) and niij) are the number of zeros and ones in the dyadic labelling 
of the state i. 

From the random /3-model all one obtains is a statement about these prob- 
abilities. This suffices to interpret most of the current experimental results 
which concern mostly velocity increments over a distance r and turbulent 
energy dissipation over a box of size r. These quantities do not involve time- 
correlations and therefore do not make a detailed test of the dynamics of the 
cascade, they only test its invariant probability measure. As discussed above 
to identify the physical mechanisms behind the structure of fully developed 
turbulence, more information is needed. If one wants, for example, the time 
correlations at a point moving with the free-stream velocity of the fluid one 
should explicitly consider models for the dynamics in state space at each 
level n. To the same invariant measure pi correspond many different pro- 
cesses. The most unstructured process corresponds to the statement that, if 
at time zero one finds the state i, then the transition probability to the state 
j at time t is proportional to pj. For the turbulent cascade the unstructured 
process does not seem to be natural because, if the lifetime of the eddies in 
the inertial range scales like -h-^ then we expect larger eddies to live longer 
than small eddies. That is, if at time t the fiuctuation 5vn{x) at the point x 
is receiving its energy through a fragmentation history leading to the state 
i then, a short time thereafter, we expect to find a different state which is 
nearby in the sense of the natural ultrametric distance in the tree. 



To characterize a stochastic process on a tree one has to solve the Chapman- 
Kolmogorov equation for the transition probabihties 

dtpizt\yO) = [ dx{W{z\xt)p{xt\yO) - W {x\zt)'p{zt\y^)} (1.7) 



with kernels W^(z|xt) that reflect the (natural) ultrametric distance in the 
tree. For kernels that depend only on the distance VF(2;|xt) = W{\z — 
x|), Ogielski and Stein^^ found the solution of Eq.(1.7). Albeverio and 
Karwowski^"^'^^ have also constructed the stochastic processes on arbitrary 
p-adic fields Q^ for the case where the jumping kernels depend only on the 
distance between p-adic balls (see also Brekke and Olson^^). However it is 
easy to see from the equation for the probability densities 

dtp{zt) = f dx{W{z\xt)p{xt) - W{x\zt)p{zt)} (1.8) 

that if W^(z|xi(:) = W{\z — x\) then the invariant density is p(z)=const. For 
the stochastic process of the turbulent cascade we require a non-constant 
invariant density as in Eq.(1.6) and the results of the authors of Refs. 13-16 
cannot be used. 
From (1.8) if follows that with 

W{z\xt) = p{z)f{\z-x\) (1.9) 

the invariant density is p(z) and, at the same time, full account is taken of the 
dependence of the transition probability on the distance between the points 
z and X in state space. In the next Section we characterize the solutions of 
the Chapman-Kolmogorov equation for kernels of the form (1.9). In Section 
3 we then show how to use these solutions to compute (or parametrize) the 
time correlations of the turbulent cascade. 



2 Random walk on a tree with asymmetric 
jumping kernels 

We rewrite Eq.(1.8) in matrix form 

|p(t) = Wp{t) (2.1) 



where W is the matrix 



(Wii 


PlCl 


Pi £2 


Pie2 


Pies 


Pies 


Pies 


Pies 






p2ei 


14^22 


P2e2 


P2e2 


P2es 


P2e3 


P2e3 


P2es 






P3e2 


P3e2 


W^ss 


Psei 


Pses 


Pses 


Pses 


Pses 






P4e2 


P4e2 


P4ei 


W,, 


P4es 


P4es 


P4es 


P4e3 






Pses 


Pses 


P5es 


pses 


W^55 


Psei 


P5e2 


P5e2 






Pees 


Pees 


Pees 


Pees 


p2ei 


M^ee 


P6e2 


P6e2 






PtCs 


p7es 


Pres 


P7es 


P7e2 


P7e2 


Wjj 


Prei 






Pses 

v . 


Pses 


Pses 


Pses 


Pse2 


Pse2 


Psei 


Wss 







Notice that the matrix has increasingly larger non-diagonal blocks of size 

2* X 2* which have the common factor e^. These blocks correspond to jumps 

to a ultrametric distance i. A matrix element in a block of size 2* x 2* at 

line k equals p^ej, e^ being the value of f{\z — x\) in Eq.(1.9) for a jump to a 

distance i. The elements Wu in the diagonal are such that the columns add 

to zero. 

As in the symmetric case studied by Ogielski and Stein^^ we find the complete 

set of eigenvectors of the matrix W. For a matrix of dimension 2", which 

describes the stochastic process at the nth level of the tree, the eigenvectors 

are: 

(i) The eigenvector 

( Pi\ 

P2 



\P2"/ 

with eigenvalue Aq = 0; 

(ii) n classes with 2""^^ (k=l ... n) eigenvectors each, where each eigenvector 
has only 2^ non-zero elements. The first 2^^"^ elements are positive and 
the others are negative. The 2*^ non-zero elements of an eigenvector have a 
common ancestor, in the tree, at the level n-k. The non-zero elements of an 
(unnormalized) eigenvector are formed by multiplying the corresponding pi 
by the sum of the p/s of the complementary group in the non-zero set of 
elements. The formation rule is easier to understand from an example. 
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Let n=3 (Fig. 2). There are then three classes of eigenvectors in the group 

(ii), typical examples of which are: 

a) 

/ PIP2 \ 
-P2P1 










V 



/ 



Ai = -(ei(pi + P2) + e2(p3 + Pa) + ealps + pe + p7 + ps)) 

Four eigenvectors of this type corresponding to the independent groups of 

two elements with a common ancestor at level 2. 

b) 

/ Pl{P3+ Pi) \ 

P2{P3+Pa) 

-P3{pl+ P2) 

-pAi.pl + P2) 




V / 

A5 = -(e2(pi + P2 + P3 + Pa) + eslps + P& + p7 + Ps)) 
Two eigenvectors of this type. 

c) 

/ Pl(P5+P6+P7 + P8) \ 

P2(P5+P6+P7 + P8) 

P3(P5+P6+P7 + P8) 

P4(P5+P6+P7 + P8) 

-P5(pl +P2 + P3 + P4) 

-P6(pl +P2 + P3 + P4) 

-p7(pl + p2 + p3 + P4) 

V-P8(pl +P2 + P3 + P4)y 

A7 = -e3(Pl + P2 + P3 + P4 + P5 + P6 + P7 + Ps) = -^3 

One eigenvector of this type. 
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The rule of formation for the eigenvalues is clear from the example above. 
Each eigenvalue is a sum of terms 

n 

^^ = -Ee.EPfc (2-2) 

j=l{i) k 

where eni) is the first e which covers all the non-zero elements of the vector 
and the sum X]fc Pk contains all the probability densities of the states reached 
by an e^-jump. 
A solution of Eq.(2.1) is an arbitrary superposition 

p{t) = Y,c,e->^^'vi (2.3) 

of the eigenvectors above. From (2.3) it is now easy to construct the general 
solution of the Chapman- Kolmogorov equation. To obtain the transition 
probability p{jt\iQ) from the state i at time zero to the state j at time t, one 
chooses the coefficients q in (2.3) in such a way that, at time zero, only the 
i-th component is non-zero and then read the value of the j-th component 
at time t. Before stating the general result we illustrate it by writing the 
transition probabilities for transitions between typical states in each group 
for the case n=3. 



p(lt|lO) = (^ + f^^ e-t{^i{pi + P2) + £2(^3 + Pa) + e3(p5 +'" +P8)} 

P1+-+P8 Pi{pi + P2) 



I Pi(p3 + P4) g-t{e2(pi +•■■ +Pa) + e3(p5 +■■■ +P8)} 

(Pi + P2)(Pl + P2 + P3 + Pa) 



^ Pl(P5+-+P8) ^- 

(Pl+-+p4)(pi+-+p8) 



^t€3 



p{2tm = -^ P^P^.-t^. + P2{p^ + Pa) ,,, 

Pi +••• +P8 Pl(Pl+P2) (Pl +P2)(Pl +P2 + P3 + P4) 

^ P2(P5 +- +Ps) ^-tXr 

(Pl+-+P4)(pi+-+P8) 



Pi +••• +P8 (Pl +P2)(Pl +P2 +P3 +P4) (Pl +••• +P4)(Pl +••• +P8) 

p(5t| 10) = ^^ ^^^^^ I; ^P'^ -e-*^^ (2.4) 

P1+-+P8 (Pl+-+P4)(pl+-+P8) 

Of course in this case pi +•• +p^ = 1 but we have kept this term to emphasize 

the rule of formation of the coefficients. The general rule for the transition 

probability p{jt\iO) between two states i and j at the level n in the tree is the 

following: 

(i) p{jt\iO) is a sum of terms, the first of which is pj (the target probability), 

and has as many terms as the number of eigenvectors that have non-zero 

elements in both the i and the j positions. 

p{jtm=p,+J2cke^'^'' (2.5) 

k 

(ii) The exponential factor in each term contains the eigenvalue of the asso- 
ciated eigenvector. 

(iii) The coefficients all contain in the numerator the target probability pj 
multiplied by the sum of the probabilities of non-zero elements of the corre- 
sponding eigenvectors in the half that does not contain j. The denominator 
is the sum of the half that contains i multiplied by the sum of all probabilities 
associated to the non-zero elements of the eigenvector, 
(iv) The sign of the coefficient is the product of the signs of the i and j entries 
in the eigenvector. 

3 Time correlations. Application to the tur- 
bulent cascade 

Once the transition probabilities p{zt\yO), solutions of the Chapman-Kolmogorov 
equation are known, the time correlations of the process are obtained from 

(x(t)x(O)) = / dydx y p{yt\xO) x p{x) (3.1) 



or 

(x(t)x(O)) =^Xjp{jt\iO)xipi (3.2) 

for a discrete state space. 

Using the results of Section 2 (Eq.(2.5)) one sees that at large times the 
time correlation at level n will be dominated by the largest non-zero eigen- 
value A2"-i = —en- Assuming that the dynamics of the turbulent cascade 
is controlled by the decay of the eddies, the largest non-zero eigenvalue will 
always be the same, associated to the mean lifetime of large eddies. However 
the asymptotic long-time correlation will be difficult to measure because of 
the small values of (x(t)x(O)) at large t. Error bars, in numerical or actual 
experiments, are likely to be larger than e"*^"* for t large. 
If, as we are proposing, the time correlations in the turbulent cascade are 
described by a stochastic process with kernels that depend on tree distances, 
a first qualitative prediction is the occurence of several exponential slopes, 
as the time increases, in the time-correlation functions. Notice that the 
existence of different time scales, as a consequence of the advection of small- 
scale eddies by large-scale motions, was already pointed out by Kolmogorov 
(see Ref.9). 

Of special interest is the slope of the short-time correlation which is controlled 
by the smallest eigenvalue. Using the dyadic expansion to label the points 
Xi in state space 

X, = /?o"««/5r(^) (3.3) 

p. = 7"oW(i_^)ni« (3,4) 

where no(i) and ni(i) are the number of zeros and ones in the dyadic expan- 
sion of i. Assuming 7 > (1 — 7) the smallest eigenvalue for the dynamics at 
level n is 

AS'^^ = -{eS")(pS")+prV4"^(pr^+pl"V4"^(pi"^+---+pi"V4"^(p^"^+---+pS6V } 

(3.5) 
If the dynamics of the turbulent cascade is associated to the decay of the 
eddies of different sizes, it is reasonable to assume that 

^"'^ = ^1 (3.6) 



Using this relation and the relations between the probability densities at 
the levels n and n-1 one obtains 

One concludes that the ratio of short-time correlations measures the differ- 
ence between the lifetimes of the structures at different length scales. If the 
dynamics of the cascade is controlled by the decay of the eddies and these 
have different lifetimes at different scales, the ultrametric stochastic model 
is an appropriate way to parametrize the dynamics and to characterize it in 
quantitative terms. Other models yield different correlation structures. 
Notice that here we are concerned with the time fluctuations of the turbulent 
cascade itself, not with the changes induced by the overall motion of the 
fluid. This means that for a fluid in motion with free-stream velocity U the 
correlations to measure, for an observable A, are 

{A{x + Tft,t)A{x,0)) 

The measure of the short-time behaviour of such quantities and the detection 
of several time scales in the time-correlations will test the usefulness of the 
turbulent cascade process proposed in this paper. Notice however that, in 
particular, the accuracy needed to detect different time scales, is a great 
experimental challenge. 

References 

[1] A. N. Kolmogorov; Dokl. Akad. Nauk SSSR 30,301; 31,338; 31,99; 32,16 
(1941) 

[2] B. B. Mandelbrot; J. Fluid Mech. 62 (1974) 331. 

[3] U. Frisch, P. Sulem and M. Nelkin; J. Fluid Mech. 87 (1978) 719. 

[4] A. N. Kolmogorov; J. Fluid Mech. 13 (1962) 82. 

[5] F. Anselmet, Y. Cagne, E. J. Hopfinger and R. A. Antonia; J. Fluid 
Mech. 140 (1984) 63. 

[6] C. Meneveau and K. R. Sreenivasan; J. Fluid Mech. 224 (1991) 429. 

10 



[7; 



[9 
[10 

[11 

[12 
[13 
[14 

[15 

[16 



U. Frisch and G. Parisi; in "Turbulence and Predictability in Geophys- 
ical Fluid Dynamics" (ed. M. Gil, R. Benzi and G. Parisi) pp. 84-88, 
North Holland, 1985 

R. Benzi, G. Paladin, G. Parisi and A. Vulpiani; J. Phys. A: Math. Gen. 
17 (1984) 3521. 

J. C. R. Hunt and J. C. Vassilicos; Proc. R. Soc. Lond. A434 (1991) 183. 

R. H. Kraichnan; Phys. Fluids 9 (1966) 1728. 

J. C. H. Fung, J. C. R. Hunt and R. J. Perkins; in "Advances in Turbu- 
lence", A. V. Johansson and P. M. Alfredsson (Eds.), vol. 3, Springer, 
1991. 

T. S. Lundgren; Phys. Fluids 25 (1982) 2193. 

A. T. Ogielski and D. L. Stein; Phys. Rev. Lett. 55 (1985) 1634. 

S. Albeverio and W. Karwowski; "Diffusion on p-adic numbers". Pro- 
ceedings of Nagoya Conference, ed. T. Hida, 1991 

S. Albeverio and W. Karwowski; "A random walk on p-adics - The 
generator and its spectrum", Bochum preprint SFB 237/no.l35 (1991) 

L. Brekke and M. Olson; "p-adic diffusion in glasses", Chicago preprint 
EFI-89-23 (1989) 



Figure captions 

Fig.l The state space at level n for a dyadic turbulent cascade 
Fig. 2 The three types of stochastic transitions associated to three different 
classes of eigenvectors 
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